Motion of nanodroplets near edges and wedges 
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Abstract 

Nanodroplets residing near wedges or edges of solid substrates exhibit a disjoining pressure 
induced dynamics. Our nanoscale hydrodynamic calculations reveal that non-volatile droplets are 
attracted or repelled from edges or wedges depending on details of the corresponding laterally 
varying disjoining pressure generated, e.g., by a possible surface coating. 
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Wetting phenomena 



DD 



and thin film dynamics 



on rough ^] 



or topographically 



•apm< 

fa 



and 



structured |5j, |6f substrates have been studied in detail both experimentally 
theoretically |4], III [h]] . To a large extent these investigations are motivated by the fact that 
generically substrates are rough, in particular those playing an important role in techno- 
logical processes such as oil recovery, coating, lubrication, and paper processing. A second 
driving force for the progressive application of these micro- and nanofluidic processes stems 
from the lab-on-a-chip concept which integrates pipes, pumps, reactors, and analyzers into 
a single device allowing for a cost efficient handling of minute amounts of liquid containing, 
e.g., DNA or proteins |l 1| . 

Most theoretical investigations of wetting of structured substrates have been concerned 
with thermal equilibrium. Only recently efforts have been made towards understanding the 
corresponding dynamics. Using various numerical techniques, time-dependent free surface 
flow over topographic features has been investigated I2I lie 
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151 ]. However, the appli- 



cability of these results to the nanoscale is impeded because the effect of intermolecular 

nn 

interactions, relevant at the nanoscale, either has not been considered 12, 13 or only in a 

nn 

rather crude way [U, H5J] . 

In contrast, here we study nanoscale fluid dynamics on topographical surface structures 
by properly taking into account the spatial variation of the long-ranged intermolecular in- 
teractions. We focus on edges and wedges as two paradigmatic geometric structures and 
provide a detailed dynamical study of their effects on nanodroplets positioned in their vicin- 
ity. We consider a partially wetting, non-volatile, and incompressible liquid film composed 
of a nanodroplet on top of a precursor wetting layer on a solid substrate. We take into 
account that the surface of the edge is possibly covered with a thin solid coating layer of a 
material different from the bulk substrate which allows one to tune the substrate potential 
[3] . Our analysis shows that the dynamics of droplets on such structures depends on subtle 
details of the substrate properties. 

For given intermolecular interactions, based on density functional theory the free energy 
of a prescribed liquid film configuration in contact with, e.g., a wedge-like substrate is a 
functional of the liquid-vapor interface shape 17[ which can be expressed in terms of the 
so-called effective interface potential uj [itI . Il8 |. For a flat substrate and film thickness h 
one has u(h) = — J h °° H.(y) dy where II is the disjoining pressure (DJP). Since the equation 
of motion can be expressed in terms of II, we determine II directly following Ref. |lfi| . 
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Assuming Lennard- Jones type pair potentials V a p{r) = M a p/r 12 — N al 3/r 6 , where M a p and 
N a p are material parameters, and a and (3 relate to liquid (/), substrate (s), or coating (c) 
particles, the DJP corresponding to chemically homogeneous substrate is given by 

II(r)=/ [fiV u {T-Tr)- Pl p t V al {v-i>)]tfT> (1) 
Jn s 

with r = (x, y, z) G M 3 and p\ and p s are the number densities of the liquid and the substrate, 
respectively Q s is the substrate volume. Equation (J2J) assumes that the vapor density is 
so low that its contribution to II can be neglected. For a non-coated edge Q s = {r £ M 3 | 
i,t/<0,z£l} (see Fig. |TJ) this yields 



,,12 . , 

n s | r — r' | Jn s |r - r' 

where AM = pfMu — pip s M s i and AiV = pfNu — pip s N s i. The first (second) term dominates 

near (far from) the substrate. On a flat and homogeneous substrate the equilibrium thickness 

of the wetting layer is given by the zero of II which corresponds to the minimum of u. For 

the interactions considered here, AM > is a necessary condition for the occurrence of an 

equilibrium wetting layer of nonzero thickness. Both integrals in Eq. (J2J) can be calculated 

analytically and one obtains the DJP as the corresponding difference of two contributions 

Il e = Ilg — rig 2 . In order to enrich the model we consider in addition the case that the 

substrate is covered by a thin coating layer of thickness d. Actual coating layers have a 

more complicated structure, in particular around edges or wedges, which depends on the 

specific combination of coating and substrate material as well as how the coating is deposited. 

Such details do influence the motion of droplets very close to the edge but corresponding 

calculations carried out by us indicate that a simple model with rather small d captures 

correctly the dynamics at lateral distances from the edge larger than d. The contribution of 

a thin coating layer to the DJP can be determined as above, assuming a van der Waals type 

interaction between the coating and the liquid particles. We do not consider an additional 

repulsive part of the liquid-coating interaction as this is shorter ranged (~ y~ 10 ) than the 

corresponding part Ilg 2 ~ y~ 9 lallSl- The contribution II" (x, y) of the coating layer on the 

horizontal upper part of the edge {(x, y, z)\x < 0, y = 0} to the DJP can also be calculated 

analytically. By symmetry, the contribution of the vertical coating layer is U^{y,x). Thus 

the DJP with the coating is 

n ce (x, y) = n e (x, y) + U u c (x, y) + U u c (y, x). (3) 
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FIG. 1: Edge configuration. A nanodroplet (full line) is placed near an edge and thus exposed to 
the action of the laterally varying DJP. Contour plots of the corresponding dimensionless DJP for 
the minus ((a), B = —1, C = 1) and plus ((b), B = —2.5, C = 1) case are shown (see Eq. Q). 
The edge surfaces are at x* = and y* = 0, respectively. 

Far from the edge, i.e., for x — > — oo, the DJP reduces to that of a coated /lat substrate: 
U cf (y) = vrAM/(45?/ 9 ) - vrAiV/(6y 3 ) + nAN'd/2y' 1 , with AN' = pfN u - pip c N cl measuring 
the interaction strength of the coating layer. We introduce dimensionless quantities (marked 

1 /6 

by *) such that lengths are measured in units of b = [2AM/ (15 |AiV|)] which for AiV > 
is the equilibrium wetting film thickness on the uncoated flat substrate. For the relation 
between b and the equilibrium wetting film thickness on the coated substrate see the insets 
in Fig. El The DJP is measured in units of the ratio a/b where a is liquid- vapor surface 
tension. Thus the dimensionless DJP U*j- = U c fb/a far from the edge has the form 

K f (y*) = c^±^ + ^). (4) 

The dimensionless amplitude C = Ab/a, where A = 7r(AM/45) _1 ^ 2 (| AN\/6) 3 ^ 2 , compares 
the strengths of the effective intermolecular forces in the uncoated case and of the surface 
tension forces. The amplitude B = nAN'd/(2 A 6 4 ), which can be positive or negative, 
measures the strength of the coating layer. Note that an analysis of the DJP, which is more 
refined than Eq. yields 5^0 even in the absence of a coating layer P, Q|; therefore 
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in the following we consider B as an independent parameter. In the second term on the 
right hand side of Eq. (|3J) the upper (lower) sign corresponds to AiV < (AiV > 0). In 
the following we shall refer to these cases as the plus and the minus case (see Fig. |2J). In 
the minus case, the DJP has a single zero for all B, while in the plus case there is no zero 
for B > —1.57 and two zeros for B < —1.57. A DJP with only one zero corresponds to 
a substrate which can undergo a second-order wetting transition, while two zeros indicate 
a possible first-order wetting transition. The absence of a zero corresponds to complete 
wetting 2]. The dimensionless form of the DJP for a coated edge is given by 



The physically possible ranges of C and B in these cases can be inferred from considering 
the macroscopic equilibrium contact angle 9 eq given by cos^eq = 1 + oj^iyl) with y$ being 
the minimum of the corresponding effective interface potential oj*f(y*) = j^Yi* c ^{y') dy' j2j. 
|cos#eq| < 1 implies B < —1.87 for the plus case so that the disjoining pressure has two 
zeros. In the minus case there is no limitation for B. 

Figures^a) and (b) show typical examples of the DJP at an edge for both the minus and 
the plus case, respectively. If a droplet is placed near the edge, it is exposed to the lateral 
gradient of the disjoining pressure resulting in a lateral force on the droplet. 

In the following we restrict our analysis to ridges translationally invariant along z, 
expecting that our arguments and conclusions carry over qualitatively to actual three- 
dimensional droplets. For a ridge, the lateral force density on the droplet is /* = 
t^t f gQ * Il* e (x*, y*) n x * dS*. dfl^ and are the dimensionless droplet surface and volume, 

d d 

respectively, and n x * is the x*-component of the unit surface normal vector pointing out- 
ward, n. In Fig. El this force density, estimated by the force density on a parabolic droplet 
crossing over smoothly to a wetting layer of thickness y$ (discussed as initial configuration 
for the Stokes dynamics below, see also Fig. [TJ, is plotted as a function of the distance w* 
(see Fig. 0J) between the three-phase contact line and the edge for different values of B. For 
the minus case there is a critical value B c ~ —10 which depends weakly on the drop size. 
For B > B c the force is always positive and increases towards the edge. Thus one expects 
the droplet to move towards the edge. However, for B < B c the force changes sign from plus 
to minus upon approaching the edge and the droplet will stop at a certain distance from the 
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edge which increases with decreasing B. In the plus case the force is always negative, i.e., 
pointing away from the edge, and increases upon approaching the edge. Therefore in this 
case one expects a droplet to move away from the edge. 

We analyze the liquid flow of this motion in terms of a two-dimensional Stokes equation. 
In dimensionless form the continuity and Stokes equation read 

V • u* = , V 2 u* = - V (p* + IT) , (6) 

KJ 

where u* = {u* x * , «*» ) is the velocity vector and p* is the hydrostatic pressure. The velocity 
and time scales are Ab/fi and p/A, respectively, with the viscosity /x. Lengths and pressure 
have been scaled with b and a/b, respectively. At the liquid- vapor interface mass conservation 
yields v* = n ■ u* for the normal interface velocity v*. At the liquid-solid interface a 
no-slip condition is applied and there is no flux into the impermeable substrate. Along 
the liquid-vapor interface the tangential stresses are zero (neglecting the viscosity of the 
vapor phase) and normal stresses are balanced by the pressure, the DJP, and the surface 
tension j^tj (in dimensional form r ■ n — ph = anh with the stress tensor r and the mean 
curvature k). There is no flux through the boundaries of the box used for the numerical 
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FIG. 2: Rescaled DJP far from the edge in the plus (a) and the minus (b) case for different B (see 
Eq. (HJ). The insets show the positions of the zeros of II* j as a function of B. In the plus case II*j 
has zeros only for B < —1.57. The full curves in the insets correspond to the equilibrium wetting 
layer thickness. 
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FIG. 3: Normalized DJP induced lateral force densities /* acting on an initial droplet configuration 
(see Fig. ^) of height a* = 15 and at distance w* (see Fig. 0} from an edge for the minus (a) and 
the plus (b) case. 



calculations, guaranteeing mass conservation. We solve these equations numerically with a 
standard biharmonic boundary integral method jscj for initial droplet shapes of the form 
y*(x*; t = 0) = yl + a*{l — [(x* — g*)/a*) 2 ]}\ x *~ 9 *\ m+1 . a* is the droplet height in the center 
and half the base width, and g* is the distance of the droplet center from the edge. In this 
study m was chosen to be 10. At t — the droplet is positioned with its three-phase contact 
line (x* = —(g* — a*), y* = y^ z*) at a distance w* = g* — a* away from the edge (see Fig. 

ID). 

Figure |3] illustrates the difference of the motion of the droplets in the minus and the plus 
case for values of B and C such that 8 cq = 90° is maintained. (B and C can be varied 
such that 9 eq remains constant. In this sense the present dynamics is complementary to 
the motion of droplets caused by chemically generated contact angle gradients 2l|.) The 
droplet moves towards the edge in the minus case and away from the edge in the plus case. 
The motion towards the edge in the minus case accelerates but comes to a stop when the 
right three-phase contact line of the drop reaches the edge. In the plus case, the drop motion 
decelerates but does not stop. 

As a second system we consider a wedge-like substrate. A wedge can be viewed as being 
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composed of a flat vertical substrate and an edge. Thus, combining their contributions yields 



it 



U* cf (x*) -lTf (y*,x* 



+ Il* e (-x*,y*) + I%*(-x*,y*). (7) 

For a non-coated wedge, within our model this is in accordance with Ref. 3|. Typical 
examples of the DJP in the wedge for the minus and the plus cases are depicted in Fig. |SJ 
For different values of B the DJP induced lateral force on a parabolic ridge is shown in 
Figs.Efa) and (b) for the minus and the plus case, respectively. The force is always positive, 
i.e., pointing away from the wedge for the minus case with its strength decreasing with 
distance. Thus, one expects a droplet to move away from the wedge. For the plus case the 
force changes sign from negative to positive near the wedge at a distance which increases 
with decreasing B. Thus one expects the droplet to move towards the wedge and then to stop 
before reaching the wedge. This is indeed what is observed in the numerical calculations of 
a a liquid ridge as shown in Fig. The droplet moves away from the wedge with decreasing 
speed in the minus case and towards the wedge with increasing speed in the plus case. 
However, in the plus case the droplet stops before it reaches the wedge. Onlv if the droplet 
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FIG. 4: Motion of a droplet (a* = 15, w* = 10) near an edge for the minus (lower left, B = —1, 
C = 1.27) and the plus (lower right, B = —2.5, C = 4.23) case. For see the insets in Fig. [21 
The upper figure shows the initial droplet shape. Shown are interface profiles at t* = 339 (dashed) 
and t* = 20893 (solid) for the minus case and at t* = 210 (dashed) and t* = 93699 (solid) for the 
plus case. The chosen values of B and C provide 9 eq = 90°. 
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FIG. 5: Wedge configuration. A nanodroplet is placed in front of the wedge and its motion under 
the action of the laterally varying DJP is followed. Contour plots of the DJP of the system for the 
minus ((a), B = —1, C = 1) and the plus ((b), B = —2.5, C = 1) case are shown. 

is driven into the wedge (e.g., by external forces) or if it touches the wedge during the 
initial relaxation process, the droplet gets trapped in the wedge and forms a configuration 
symmetric with respect to the bisector of the wedge. 

In summary we have shown how topographic substrate features generate motion of non- 
volatile nanodroplets residing in their vicinity. This motion depends on the details of the 
interplay between the liquid-liquid and the liquid-substrate interactions. On substrates for 
which in the planar case the DJP approaches zero from above at large distances (plus case), 
droplets move towards edges and away from wedges, and vice versa in the opposite case 
minus case). This resembles similar phenomena occurring near chemical heterogenities 

Taking 1 nm and 0.02 N/m as typical values of b and a, respectively, and for the values 
of C and B chosen in Figs. H] and our model calculations predict that near an edge the 
DJP induced lateral force on a liquid ridge 30 nm long, 15 nm high, and 1 nm thick is of 
the order of 10 -13 N. For this droplet size the gravitational force due to a possible tilting of 
the substrate is about eight orders of magnitude smaller. Only for micron-sized drops the 
gravitational body force becomes competitive. The time scale fi/A for the motion is roughly 
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FIG. 6: Normalized DJP induced lateral force densities /* acting on an initial droplet configuration 
of height a* = 15 and at distance w* (see Fig. [7J) from a wedge for the minus (a) and the plus (b) 
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FIG. 7: Motion of a droplet (a* = 15, w* = 10) near a wedge for the minus (lower left, B = — 1, 
C = 1.27) and the plus (lower right, B = —2.5, C = 4.23) case. For see the insets in Fig. |2I 
The upper graph shows the initial droplet shape. Shown are interface profiles at t* = 450 (dashed) 
and t* = 99750 (solid) for the minus and t* = 160 (dashed) and t* = 14683 (solid) for the plus 



case. The chosen values of B and C provide 
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[lCT 8 /i/ (Pa s)] s so that the average velocities for the cases studied are ca. {lCT 5 /[/i/(Pa s)]} 
m/s. Taking // between 0.1 Pa s and 100 Pa s (different Polydimethylsiloxanes at ambient 
temperature) the velocity ranges from 0.1 mm/s to 0.1 //m/s, which is comparable with the 
velocities of droplets of the same kind of liquid but exposed to and driven by a chemical 
step generated by a contact angle contrast of a few degrees |2l|. Our calculations indicate 
that the initial relaxation of the droplet shape from the prepared one, i.e., spreading of 
the droplet radius, occurs faster than the subsequent motion of the center of gravity of 
the droplets. Since the spreading rate of the monolayer precursor film outside the droplet 
tends to be larger than that of the apparent droplet radius our calculations lead to the 
expectation that the velocity for the precursor spreading is larger than the motion of the 
droplets. 

Our findings have a bearing on the distribution of droplets after condensating liquid onto 
nano-sculptured substrates and on breath figures, in addition to possible implications for 
the design of open nanofluidic systems in which the activating force, partially or totally, is 
provided by the system itself. 
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